Local edge modes in doped cuprates with checkerboard polaronic heterogeneity 
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We study a periodic polaronic system, which exhibits a nanoscale superlattice structure, as a 
model for hole-doped cuprates with checkerboard-like heterogeneity, as has been observed recently 
by scanning tunneling microscopy (STM). Within this model, the electronic and phononic excitations 
are investigated by applying an unrestricted Hartree-Fock and a random phase approximation (RPA) 
to a multiband Peierls-Hubbard Hamiltonian in two dimensions. 



Recent scanning tunneling microscopy (STM) studies 
of a Cu oxychlorideQ, Q > as wen as some other doped 
cupratesjj, 0, M, M, that exhibit high-temperature su- 
perconducting behavior, has revealed that doped holes 
can drive mesoscopic domains of novel nanoscale elec- 
tronic modulation at low doping [x ~ jg), that is, 
"checkerboards" , rather than a simple stripe modulation. 
The checkerboard domains, according to analysis of the 
experimental data|l|,|2j, consist of (4 x 4)-unit-cell pla- 
quettes; the plaquette size is essentially insensitive to the 
doping level, presumably locking to this commensurate 
periodicity with discommensuration defects. 

However, it is not clear where doped holes are located. 
If the holes lie along the grid of a checkerboard, the linear 
hole concentration on the grid would be extremely dilute 
for underdoped cases. It is unlikely, taking into account 
that the system is an insulator at such a low doping level, 
that holes compose a 4-site-periodic checkerboard and are 
diluted throughout the whole plane like a metal. Thus 
it seems more likely that each doped hole is localized, 
and forms a "polaron" with associated lattice and spin 
distortions. 

Here, therefore, we study a periodic polaronic 
structure as a model for hole-doped cuprates with 
checkerboard-like superlattice heterogeneity. Our aim is 
to determine the spectroscopic signatures of such checker- 
boards, and contrast them with signatures of stripe seg- 
ments. We expect that the results of this work will aid 
in the analysis and interpretation of recent experimen- 
tal data of angle-resolved photoemission spectroscopy 
(ARPES), nuclear quadrupole resonance (NQR), elec- 
tron paramagnetic resonance (EPR) , an d the dynamical 
structure factor S(q, lo) ^,%,M,m^,^Ek 

There have been several attempts to describe cuprate 
superconductors as a system with stripe modulation. 
We|l5j previously investigated the phonon mode struc- 
ture in the presence of stripes via Hartree-Fock and RPA 
analysis in a three-band Peierls-Hubbard model, and re- 
ported that, owing to the stripe inhomogeneity, local 
edge (interface) modes of oxygen lattice vibrations are 
generated with characteristic energies around 15-45 and 
70 meV; these energies are similar to those where the 
electronic band dispersion has a "kink" in ARPES 8] 
and where anomalous inelastic neutron scattering has 
been identinedpjj |l7J. These modes are also accompa- 



nied by correlated electronic excitations in both mag- 
netic and charge channels [l5j. The dispersive phonon 
branch at around 15-40 meV is due to oxygen vibrations 
along a vertical stripe and is only present for the verti- 
cal/horizontal stripe orientation. The 70 meV mode is 
from transverse oxygen vibrations and is present for all 
stripe orientations. 

In this study, we adopt a system with periodic po- 
larons forming a 4 x 4 superlattice for a -^-hole-doping 
case, and investigate the local electronic and phonon 
modes to complement and discriminate from the stripe 
cases. Although the hole-doping level here is extremely 
low {x = jg), we compare the results with those of the 
stripe cases with a different doping level {x = 4), be- 
cause it is numerically difficult to calculate such a large 
system. However, this comparison is justified, because 
the energy of the phonon edge modes in the stripe cases, 
being a local property, is essentially insensitive to the 
doping level [15j. 

To model the Cu02 planes of doped cuprates, we 
use a two-dimensional (2D) three-band extended Peierls- 
Hubbard Hamiltonian, which includes both electron- 
electron and electron-phonon interactions |l8l llflj : 

H Q = Y^ hd (Uij ) (4, c ja +H.c.) + ^2e 2 (u y )cJ CT c la 

<ij>(7 i,a 
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Here, c ia creates a hole with spin a on site i; each site 
has one orbital (d x 2_ y 2 on Cu, or O p x or p y on O). The 
Cu (O) site electronic energy is e<j (e p ). U%j represents 
the on-site Cu (O) Coulomb repulsion, Ud (U p ), or the 
intersite one, U p d, and the summation Y^ is taken, except 
for the case of (i, a) = (J, a'). The electron-lattice inter- 
action causes modification of the Cu-0 hopping strength 
through oxygen displacement Uif. t p d{uij) = t p d ± auij, 
where +(— ) applies if the Cu-0 bond shrinks (stretches) 
for a positive Uij\ it also affects the Cu on-site energies 
ed{uij) = £d + ftJ2j (^ zU ij)t where the sum runs over the 
four neighboring ions. Other oxygen modes (buck- 
ling, bending, etc) are assumed to couple to electron 
charge more weakly and are neglected for simplicity, but 
can be included as necessary within the same a ppr oach. 
We use the following set of model parameters 



pproacn. 



e p -e d = 4.4 eV, U d = 11 eV, U p = 3.3 eV, U pd = 1.1 eV, 
and K = 38.7 eV/A 2 , a = 5.2 eV/A, /3 = 1.2 eV/A, with 
ip f z = 1.1 cV. To approximately solve the model, we use 
the unrestricted Hartree-Fock approximation combined 
with an inhomogeneous generalized RPA for linear fluc- 
tuations of lattice, spin or charge 18J in a supcrccll of size 
N x x N y with periodic boundary conditions. 

The output of the calculation here is the inhomo- 
geneous Hartree-Fock ground state and the linearized 
fluctuation eigenfrequencies and eigenvectors. From the 
phonon eigenmodes, we calculate the corresponding neu- 
tron scattering cross section: 



S(k,uj)= f dt, 



IV 
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-ikR,(0) ikR, 



(t)\ 
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where R;(£) = R; + d; + u/(i) is the position of the Z-th 
oxygen atom expressed in terms of the location of the unit 
cell origin R°, the position within the unit cell d;, and 
the time-dependent vibrational component u/(t). For 
phonon modes with u; (t) oriented along the correspond- 
ing metal-oxygen bonds, on the O x sublattice d; = f x 
and u; = xix, and on the O y sublattice d; = |y and 
u/ = yiy. The scalar displacements can now be expressed 
in terms of the normal modes z n as xi(t) = J2 n a xi,nZ n {t) 
andyi(t) = X) n a 2/i>™ z «W- performing first- order expan- 
sion in the oxygen displacements, we obtain 



S(k,w) = X)[ ; 



2 \^ x 2 , .2i y |2 



x l u k,n 



k y\ a kn\ 



(3) 



+k x k y (e^ k *- k y^ 2 al n a\ n + c.c.) 



x- [(1 + n B )5(uj -w n ) +n B 8(uj +w n )l- 
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is the thermal population of the phonon mode n. This 
is a generalization of the usual neutron scattering inten- 
sity expression 20] for the case of phonons with a larger 
real space unit cell. We plot 5(k, cj)/|k| 2 for k-directions 
sampling longitudinal modes, consistent with the com- 
mon experimental convention. 

To investigate the neutral electronic excitations, we 
calculate the spectral function: 



Y,\^o\0(k)\^ n )\ 2 S(uj-E n -E a ), 



(4) 



where \^q) (l^n)) is t ne Hartree-Fock ground (excited) 
state whose energy is represented by E (E n ), and O(k) 
is an operator, e.g. spin S(k) or charge n(k), summed 
over Cu- and O-sites: 



0(k)=^oW(k) 



-lkdj 



A = Cu, 

o I ,o„ 



for A = Cu 
§£ for A = O x .(5) 
%y for A = O y 



The effect of an infinitesimal external field corresponding 
to the excitation ^ n can be represented by the change of 
an observable (O) in the state * = * + V^n (\v\ < !) : 



(O) Si (O)o + S(0) n , 

S(0) n CX (* |O|*n>, 



(6) 

(7) 



where (O)o is the expectation value with respect to the 
ground state. 

We first show the ground state configuration given by 
the Hartree-Fock calculation (Fig. Q^a)). The ground 
state has one polaron per unit cell of the polaronic 4x4 
superlattice. In the case of two holes per unit supercell, 
the extra hole is preferentially located at the central site 
of the one-hole-polaron supercell. Because of the symme- 
try, the central sites are likely to be important in under- 
standing both the phonon mode and electronic excitation 
even in the one-hole case. 
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FIG. 1: (a) The ground-state configuration given by mean- 
field Hartree-Fock calculation. Black and red circles repre- 
sent Cu and O sites, respectively. The radius of the circle is 
proportional to the corresponding site hole density. Arrows 
centered on circles show the magnitude and direction of spin. 
The green lines originating from O sites indicate the mag- 
nitude and direction of equilibrium O displacements in the 
presence of a polaron. (b) The phonon spectra calculated by 
means of RPA for the doped polaronic case (left) and the un- 
doped AF case (right); here the intensity is normalized by the 
maximum value, and magnified 4-fold below 74 meV. In the 
polaronic case (left), there is an extra branch (indicated by 
the red arrow), which is from a local vibration mode around 
a polaron, near 68 meV. 



The phonon spectra are then obtained by real-space 
RPA (Fig. IHb)). In these spectra, a dispersive branch 
occurs at 77-92 meV in both the doped and undoped 
cases. This branch is essentially insensitive to doping. 
However, doping causes an extra dispersionless branch 
near 68 meV. This branch corresponds to a local mode 
- the "shape mode" — due to the heterogeneity of the 
system with sharp interfaces, i.e., the local vibrations of 
the oxygens around a polaron. These results are very 
similar to those found in the stripe case 15] . However, 
note again that, in the vertical stripe case, there are two 
branches of local phonons at the stripe that correspond to 
oxygen oscillations with directions parallel and perpen- 
dicular to the stripe, while there is only one local phonon 
branch in the polaronic or diagonal stripe cases. 
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FIG. 2: The excitation spectra for charge (red), longitudinal 
spin (blue) and transverse spin (green) channels given by real- 
space RPA. The spectra for charge and longitudinal spin are 
scaled on the right axis, and that for transverse spin is scaled 
on the left axis, (a) and (b) are for the polaronic case; (c) 
and (d) for the undoped case. 



89 meV is shown in Fig. [3] (d). 
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Next, we calculate the electronic excitations within 
the same real-space RPA. The transverse spin mode has 
a non-zero-gap excitation (~ 8 meV), corresponding to 
AF spin waves, at several k-points: (tt, tt), (f,7r), (0, tt), 
etc. This is because the polaronic configuration cannot 
be expressed by a single k. The excitations at (tt, tt) in 
this channel has not only the ~ 8 meV peak but also 
peaks at ~ 51 meV and ~ 57 meV (Fig. |H(a)) for sub- 
spin waves (i.e., higher energy transverse spin excitations 
generated from low-energy spin waves, as shown in Fig.0 
(a)). These peaks also exist at other k-points. The ex- 
citation modes corresponding to the peak at ~ 51 meV 
are again similar to those found for vertical stripes; how- 
ever, the homogeneous system does not have this mode 
at (tt, tt). This mode occurs by the folding back of the 
spin wave branch. It depends on the symmetry of the 
superlattice; e.g. the two-hole system does not exhibit 
the sub-spin waves at (it, tt). 

Furthermore, the same kind of excitation appears in 
the undoped case (Fig. 12(d)), at other than (tt, tt). How- 
ever, this excitation has a single peak (unlike the polaron 
case) whose energy depends on the momentum. This is 
not surprising because this is the branch of the spin wave 
excitation originating at (tt, tt). 

Finally, the charge mode at (tt, tt) has a well-defined 
peak near 92 meV (Fig. El (a)), coupled with the corre- 
sponding phonon excitation. This mode shows the typi- 
cal (tt, tt) charge oscillation at the Cu sites, as in the case 
of undoped AF. The longitudinal spin mode at (~, tt) has 
well-defined peaks at ~ 77 meV and ~ 81 meV (Fig. 
(b)). These modes show the oscillation of S z at the oxy- 
gen sites, as in the case of undoped AF. However, unlike 
the undoped case, a weak peak also exists at ~ 89 meV 
at (tt, tt) in the charge and the longitudinal spin channels. 
One of the excited states corresponding to this energy at 



FIG. 3: The excited-state spin eigenvectors for (a) E ~ 51 
meV, (b) E ~ 68 meV, (c) E ~ 77 meV and (d) E ~ 89 
meV. The arrows in (a) represent S x , while those in (b)-(d) 
represent S z . 



At the characteristic energy of 68 meV, where the lo- 
cal phonon mode exists, the charge channel exhibits a 
major difference between the polaronic and the undoped 
AF cases (Fig. (b)). The charge mode of this en- 
ergy yields well-defined peaks in the excitation spectra 
at (0, tt), (0,5), (f,7r) and the symmetric equivalents 
of those points, and the intensity is the highest at (§,7r) 
— and also at the symmetrically equivalent points — in k- 
space. From the analysis of the excited state correspond- 
ing to this mode, we find that the charge on the polaron 
undergoes local oscillations of its amplitude, coupled with 
the local phonon mode (the so-called half-breathing mode 

EE3). 

In summary, we have studied a periodic 4x4- 
superlattice polaronic ground state as a model of the 
nanoscale heterogeneous structure recently observed in 
some hole-doped oxychlorides and other hole-doped 
cuprates. We investigated the characteristic electronic 
and lattice excitations relative to that ground state. The 
inhomogeneity of the system causes excitations distinctly 
different from those in the homogeneous case. The major 
new excitations occur in the charge channel (~ 70 meV 
mode), and are coupled with local half-breathing phonon 
modes, as well as in all the stripe cases. However, the 
oscillation patterns of this excitation are, compared with 
those directed along the stripe in the vertical stripe case, 
localized at the polaron site. Assuming the system pos- 
sesses a vertical stripe (or stripe segment) structure, we 
expect to find two kinds of excitations at around 15-45 
meV and 70 meV. On the other hand, if the polaronic 



superlattice structure exists, then there are only excita- 
tions at around 70 meV. A diagonal stripe case is similar 
to the polaronic case for charge excitations; however, in 
the transverse spin channel, there is an additional exci- 
tation at ~ 50 meV in the 4x4 polaronic case. The 
recent analysis of ARPES data supports the presence of 
these low-frequency local (edge) mode excitations, which 
are characteristics of polaronic inhomogeneity existing as 
either checkerboard or stripe (segment) domains. As we 
predict, these assignments are correlated with inelastic 
neutron scattering (for spin and charge) (16L |X7T | showing 
the existence of the modes at around 40 and 70 meV. 
NQR[H El IH and EPRQ data are likewise consis- 
tent with the assignment. The correlation among spin, 
charge and lattice signatures of these local (interface) 
modes will be essential for their definitive assignments 
and awaits further experimental study. 

Finally, we emphasize that similar signatures of 
nanoscale electron/lattice inhomogeneity can be ex- 
pected below the polaron formation temperature in many 
doped transition metal oxides (bismuthates, manganites, 
nickelates, etc.) and related strongly correlated elec- 
tronic materials. It will be important to compare these 
cases and identify which signatures correlate with the ap- 
pearance of superconductivity. 
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